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GAMMA-RAY BURST HIGH TIME-RESOLUTION SPECTRAL OBSERVATIONS 
M‘u)E WITH THE SOLAR MAXIMUM MISSION 

U. D. 0*0*1, J. P. Norr.'s*, T. L. Clin*, B. R. D*nnl*, K. J. Frost 
NASA/Goddard Spac* Flight C«nt*r 
Gr*«ob«lc, MD 20771, USA 


Result* from th* KONUS experiments on Venera 11 and 12 have 
significantly enriched yray burst phenomenology with spectral and 
temporal details for a large number of events (Matets at al. 1981). 
Their low energy spectral features in the 30 to 70 k*V range have been 
interpreted as due to cyclotron absorption and emission phenomena. He 
report high time-resolution (128 ms) spectral observations of Y-ray 
bursts with Solar Maximum Mission (SMM). We confirm th* existence of 
deficiencies of low energy photons using extrapolated mxponentlal- 
contlnuua fits from above 100 keV. Significant spectral variations, 
primarily at the lower energies, are also seen by SMM on time 
scale* • .25 s. Details of the SMM Hard X-Ray Burst Spectrometer 
(UXRBS ) along with a previous analysis of one of the events discussed 
here have been given by Dennis et al. (1982). The event April 19, 1980 
is reanalyzed and presented together with the events of April 21, 1980 
and March 1, 1981. 

Figure 1 shows the tlme-hlstorles of these events for three energy 
bands covering the 25 to 500 keV range. April 19 and March 1 were 
relative'y short events, • 4 s and ■ 2 s, respectively, compared to 
April 21, " 20 s. 

A Monte Carlo approach was utilized to simulate the HXRBS response 
to candidate spectral models, and thereby deconvolve the detector 
energy-loss spectra. Models representing thin thermal bremsst rah lung 
(TTB), thin thermal synchrotron (TTS) with and without absorption 
features, and power-law with exponential cutoff (PLE) were used to 
generate incident photon spectra. Best fits (least X' v ) were obtained 
for the whole energy range by including absorption features at low 
energies (TTS and TTB). To study the evolution of the spectra, the 
events were subdivided into intervals corresponding to fine structures 
in the temporal profiles such as rising phase, decay phase, and valleys 
between pulses. This division is illustrated in Figure 1. 

# 

Figure 2 gives the photon energy spectra for the three events for 
the various intervals. The parameters for the best fits, adopting the 
TTS model, are indicated for various Intervals of the events in Figure 
2. For the two April events, Che best fits are obtained by including 
absorption lines at low energy except for the valley Intervals. The 
weak short event of March 1, 1981 exhibited a smooth continuum without 


* 

^Physics Department, University of Maryland, College Park, M) 20742 
From a Ph. D. dissertation performed at the University of Maryland. 
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evidence of spectral 
faaeuraa. Tht power-law 
portion of tha PLE model 
provided aeeaptabla flea to 
tha low-energy ehannala 
without tha inclusion of 
absorption faaturas. Tha 
flattening in tha apactra 
evidenced by tha PLE aodel 
indicates, Indeed, a 
daprasslon or "shelf" at 
low energies - 

Table 1 gives tha 
values of the pa roasters 
for the TTS, TTB, and PLE 
aodela for intervals 
indicated In Figure 1. Tha 
paraaeters show significant 
variations during tha 
course of the events. The 
values of kT in tha TTB 
aodel have a vide dynamic 
range of 100 keV to 2 MeV. 
If one assuaes a constant 
value for the magnetic 
field strength of ~ 10 l G, 
then the critical energies 
deterained froa the TTS 
aodel Indicate a narrower 
range of temperatures, 180 
keV to 570 keV. For both 
April events, Share at si. 
(1982) have fitted a power 
law to PHA data at higher 
energies up to 10 MeV . 

Froa this analysis and 
recent results from the 
SIGNE experiment (Barat 
1983), It Is apparent that 
fast spectral evolution In 
Y-ray bursts aay create the 
appearance of spectral 
features _or^ camouflage 
actual features, and 
distort the contlnuua shape 
If spectra are Integrated 
for Intervals £ 1 s. 


a 
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TABLE 1 

* Thin Thermal Bressatrahlung w/ w/o Absorption Line 


AP1980 



C 1 

kT 

C 2 

Energy 

Width 

X 2 /v 

RSI 

1. A1±0. 03 

709113 

32.611.6 

43.511.6 

50.0114.1 

10.8/10 

DC1 

2.1910.13 

604133 

36.215.0 

55.14.7 

50.0117.6 

8.6/10 

VL1 

0.6310.57 

120011177 




7.8/13 

RS2 

3.1910.05 

781111 

36.011.0 

48.511.0 

33.012.2 

12.0/10 

DC2 

2. 50t0. 14 

496126 

27.313.8 

36.915.3 

50.0140.1 

14.7/10 

VL2 

0.7910.38 

9471347 

AP2180 



9.7/13 

RSI 

9.0610.23 

10312 

37.611.3 

450210.9 

50.0110.4 

20.4/10 

OC1 

4.74±0.19 

12814 

10.412.8 

41.819.3 

33.0126.9 

8.2/10 

VL1 

1 .03±0.35 

135135 




5.0/10 

RS2 

4.3210.13 

16614 

39.611.7 

46.211.4 

50.0112.0 

27.3/10 

DC2 

4.5410.48 

148110 

36.614.8 

46.512.8 

50.0129.0 

6.2/10 

VL2 

2.1810.64 

119120 

32.6115.1 

51.514.8 

39.0131.1 

3.3/10 




MR0181 




RS 

0.1010.02 

19591250 




14.5/13 

DC 

0.09110.035 

15141586 




11.4/13 


Page No. 4 


ORIGINAL PAGE IS 
OF POOR QUALITY 


Paper Code No. XG2-2 


# Thermal Synchro cron w / v/o Absorption Line 


AP1980 



C 1 

c c 

C2 

Energy 

Width 

X 2 /v 

RSI 

3.08*0.17 

10.3*0.4 

32.2*5.1 

53.7*5.6 

50.0*20.4 

9.2/10 

DC1 

5.33*0.31 

8. 4*0. 3 

34.1*5.1 

69.8.46.3 

50.0*13.3 

7.1/10 

VL1 

2.55*1.60 

9. 9*4. 7 




7.3/13 

RS2 

6.23*0.26 

11.5*0.4 

26.8*3.4 

51.4*4.3 

-33.0*10.1 

13.3/10 

DC2 

7.00*0.41 

6. 4*0. 3 

28.2*4.3 

41.5*5.5 

50.0*44.2 

12.9/10 

VL2 

2.77*1.10 

9. 2*3.0 

AP2180 



8.4/13 

RSI 

52.2*10.9 

1. 4*0.1 

35.8*4.6 

44.5*1.9 

50.0*26.9 

22.9/10 

DC1 

38.0*1.5 

1. 5*0.0 

17.8*2.4 

42.6*4.5 

33.0*12.4 

7.3/10 

VL1 

4.88*0.37 

1. 9*0.1 




5.0/13 

RS2 

18.6*0.6 

2. 3*0.0 

41.9*1.7 

47.3*1.3 

50.0*9.9 

31.2/10 

DC2 

4.04*0.91 

4. 4*0. 6 




12.7/13 


14.8*1.6 

2.4*0. 1 

32.4*5.4 

47.6*3.5 

50.0*32.4 

5.6/10 

VL2 

2.49*0.69 

3.340.6 




12.7/13 


13.0*4.4 

1.640.2 

34.7*12.6 

50.7*3.8 

40.9*26.7 

3.4/10 




MR0181 




RS 

0.56*0.16 

11.9*2.6 




11.1/13 

DC 

0.30*0.06 

14.5*2.4 




10.2/13 


RSI 

C 1 

1.04*0.05 

Power-law with Exponential Cutoff 
AP1980 

° E cut E char 

0.57*0.01 196*18 146*37 

X 2 /v 

4.5/11 

DCl 

3.01*1.14 

0.76*0.09 

250*90 

100*120 

5.6/11 

VL1 

2.40*0.58 

0.78*0.05 

338*163 

19*130 

7.7/11 

RS2 

1.02*0.04 

0.37*0.01 

127*9 

387*91 

16.5/11 

DC2 

0.89*0.93 

0.46*0.24 

132*37 

171*77 

9.7/11 

VL2 

4.71*8.95 

0.89*0.48 

118*111 

241*184 

8.1/11 

RSI 

8.18*2.76 

1.10*0.10 

AP2180 

130*15 

68*15 

15.4/11 

DCl 

15.2*4.8 

1.22*0.09 

73*28 

102*20 

8.8/11 

VL1 

20.1*11.5 

1.63*0.02 

58*246 

112*54 

5.0/11 

RS2 

3.46*1.50 

0.95*0.10 

193*29 

63*26 

15.2/11 

DC2 

6.19*4.32 

1.08*0.17 

129*40 

113*40 

5.2/11 

VL2 

5.46*4.48 

1.27*0.20 

124*49 

99*44 

8.0/13 
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DIFFUSE GALACilC EMISSION OF HIGH ENERGY GAMMA RAYS 

D.A. Kniffen, C.E. Flchtel, and R.C. Hartman 
NASA/Goddard Space Flight Center 
Greenbelt, Maryland 20771 U.S.A. 


1 . Introduction 

Assuming cosmic rays pervade the Galaxy, they necessarily produce 
high energy y~rays as they Interact with the Interstellar matter and 
photons. Substantial work has already been performed on the calculation 
of the source functions for these various y radiations and the intensity 
to be expected in the vicinity of the solar system. (For a general 
review see Chapter 5 of Flchtel and Trombka, 1981.) However, several 
recent devlopments make a reexamination and extension of this work 
worthwhile. These include the recently published detailed results of 
high energy galactic y-radlation obtained with the COS-B satellite 
(Mayer -Has selwander et al«, 1982) and the new COS-B background estimate 
(Strong, 1982), further evaluations of the 21 cm radiation in the galaxy 
and, hence, the atomic hydrogen density distribution, considerations 
related to the molecular hydrogen density normalisation, and improved j 

theoretical calculations on the nucleon-nucleon source function. j 

5 

2 . Gamma Ray Production and Intensity Calculations f 

$ 

The galactic matter distribution and content by species are ! 

discussed in earlier work (e.g., Kniffen, Flchtel, and Thompson, 1977), * 

and only more recent consideration will be mentioned here. The density j 

of neutral atomic hydrogen as revealed by the 21 cm emission remains 
somewhat uncertain in the inner galactic regions because of uncertainty 
in the absorption correction. Recent work (e.g. Dickey et al., 1982; j 

Thaddeus, 1982) suggests that the absorption had previously been J 

somewhat underestimated. Thus, in this work, the atomic hydrogen 
density distribution of Gordon and Burton (1976) as a function of radius 
from the galactic center was used, but modified so that the atomic } 

hydrogen density in the Innermost region was increased by a factor of ] 

1.5, and the closer densities were Increased less in accordance with the | 

amount of intervening matter. Recent work by Lockman (1982) has, shown 1 

that the scale height is about 1 1/2 times larger then previously j 

believed because a relatively faint component of HI had been overlooked. 1 

Hence, the scale height used in this work is 0.18 kpcs for a galactic j 

radius, Rq 8 i, less than that of the sun (10.0 kpcs) and [0.18 + 0-023 x 
(Rg a l - 10.0)] kpcs for Rg^ > 10.0, with the incease beyond the solar 
circle being based on the work of Baker and Burton (1975). It is now 
believed that the scale height in the outer galaxy increases more 
rapidly than this (Kulkarni, Blitz, and Heiles, 1982), but the surface 
density is still believed to be similar to earlier estimates. The 
density distribution of molecular hydrogen is measured less directly 
than that of atomic hydrogen, and hence the absolute value Of the 
density is less certain. In this work the galactic distribution of 
Gordon and Burton (1976) was used, but density normalization was treated 
as an adjustlble parameter. 

Although the translation of the 21 cm observations into a galactic 
spatial distribution is difficult, on a broad scale the density profile 
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1* reasonably veil known, even thougjh details of arm structure are not 
always agreed on by all workers in the field* Until recently, it had 
not been clear whether molecular clouds were associated with spiral 
structure. However, now on the basis of a high saapl* survey and obser- 
vations In both the first and second quadrants of the galactic plane, 
Cohen et al. (1980) have reported the existence of the isolecular coun- 
terparts of the five classical 21 ca spiral eras segments In these 
quadrants* The specific spiral pattern that will be used here Is that 
of Georgelin and Georgelln (1976). A five hundred parsec width Is 
adopted for the arms. The excess of material In the aras Is taken to be 
the smaller of 2SZ of the total local matter or twice the local average 
density of matter not In the aras, based on recent considerations 
(Lockaan, 1982, and Kulkarnl, Blits, and Kelles, 1982). 

For the photon distributions, Kniffen and Fichtel (1981), using 
results of Bolssd et al. (1981) on the Infrared volume emlssivity and a 
model of Bahcall and Sonelra (1980) for the starlight distribution, 
obtained photon densities and, hence, source function for the Compton 
emission as a function of position in the galaxy. 

With regard to the cosmic ray distribution in the galaxy, it will be 
assumed that the nucleonic cosmic ray composition and energy spectrum 
remain unchanged throughout the galaxy and that the electron spectrum 
changes only in a second order manner as the density changes. (See 
particularly Kniffen and Fichtel, s 1981, and Fichtel et al., 1976) For 
the reasons described In the above works, the cosmic ray density In the 
plane will be assumed to be proportional to the matter dcr.3ity on the 
scale of arms and, perpendicular to the plane, to have a gaussian dis- 
tribution with a scale height of 0.6 kpc . This latter number is based 
on the radio continuum measurements of Cane (1977) and the assumption 
that the galactic magnetic field energy density and the cosmic ray 
energy density have the same scale height. 

The production functions and intensity calculations are essentially 
the same as earlier work with the exception of the source function for 
cosmic ray nucleon matter interactions. Badhwar and Stephens (1977), 
Stephens and Badhwar (1981), and Morris (1982) have used the substantial 
recent high energy physics experimental work to estimate the y-ray pro- 
duction energy spectrum for cosmic rays interacting with interstellar 
matter. The spectral shapes calculated by Badhwar and Stephens (1977) 
and Morris (1982) are similar. Both are, however, significantly diffe- 
rent from the earlier pioneering work based more heavily on theoretical 
models when the more extensive experimental data were not available. 
The Morris (1982) work was used here. 

3 . Results 

The y-ray intensity predicted by this work is compared to the COS-B 
longitude distributions in Figure 1. It should be mentioned that COS-B 
scientific group has recently reanalyted the instrumental y-ray 
background (Strong, 1982). Whereas the background intensities for the 
(150-300) MeV and the (300-5000) MeV energy intervals are essentially 
unchanged from the paper of Mayer-Hasselwander et al. (1982), the (70- 
150) MeV background intensity is now estimated to be 5.6 * 10” photons 
cm“*a“* ster , a factor of -1.75 higher. This change has been incor- 
porated in Figure 1 by introducting a new "zero base". The revision 
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MLMTIC LONITUOC i 

Figure 1: Gamma ray intensity as 
a function of longitude averaged 
over the latitude range -10° < b 
< 10* from 70 feV - 150 MeV, 150 
MeV - 300 MeV, and 300 MeV - 
5000 MeV from the COS-B data 
(Mayer-Hasselwander et al., 
1982) compared to the model 
discussed heiC shown by the 
solid line. The dashed line in 
the 70-150 MeV graph represents 
the new "0.0" line based on the 
revised background intensity for 
the COS-B data in this energy 
interval discussed in the text. 

Figure 2: Energy spectrum of the 
galactic y radiation for a 
region near the galactic center. 
The calculated spectra are based 
on the work of Kniffen and 
Fichtel (1981) modified to 
include the recent nucleon- 
nucleon calculations of Morris 
(1982). The COS-B data are 
those of Mayer-Hasselwander et 
al. (1982), and the SAS-2 data 
are those of Hartman et al . 
(1979) . 

means that not only are the Intensities for the center and the sources 
(195,5), PSR 0531+21, PSR 0833-45 measured with the SAS-2 and COS-B 
instruments in good agreement with each other, but also the general 
galactic diffuse intensity. It was noted earlier in this article that 
molecular hydrogen density normalisation was left as an adjustable para- 
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meter. There did need to be an adjustment for the loner ring by a factor 
of about 0.7, but no adjustment was required beyond this ring. Regarding 
Figure 1 and the energy spectrum shown in Figure 2, considering the 
uncertainty In the point source contribution and the mass distribution, 
the agreement between the data and the predicted curves seems reasonably 
good especially when the known discrete sources are taken into account. 

The latitude distributions resulting from the model have been calcu** 
lated taking Into account the COS-B instrument response. The fit to the 
galactic center observations is remarkably good at all energies. This 
gives confidence that not only is the total observed emission well pre- 
dicted by the model, but also that the radial distribution of the emis- 
sion calculated from the model oust also be approximately correct. The 
anticenter result is reasonable considering the uncertainty of the y-r*y 
measurements and the gas distribution in that general direction. At 1*70 
to 90, the overall intensity prediction is about right, but there is some 
indication in all energy Intervals that the average emission may be 
slightly more distant than the model prediction. 

It should also be mentioned that in the concept being presented here 
the arms on the far side of the galaxy make an important contribution for 
small (|b| £ 0.4®) galactic latitudes. With future high resolution y~ ray 
measurements, these back side arms should appear as a narrow ridge super- 
imposed on the broader ridge of the near side arms . It may be possible 
to identify very large far side molecular clouds if the majority of the 
molecular hydrogen is in large clouds . 
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ELECTRON PROPAGATION IN THE LEAKY BOX MODEL 
WITH A TRUNCATED PATHLENGTB DISTRIBUTION 

B. 6* Mauger* and J. ?. Ones 

NASA/ Goddard Space Flight Cantor, Groonbolt, MD 20771, U.S.A* 
1. INTRODUCTION 


Early studios (Oraos and Proior, 1978; Protneroe, Ormes, and 
Comstock, 1981) of the choaleal coaposltlon of eosaie ray nuclei with 
3 < 2 < 28 suggostod that observed primary cosmic ray spoctra have boon 
stopponad by a faetor of E”’** _ 1 from thoir sourco spoctra, suggostlng 
source sspoctra that vary as E 2 * s “ * 1 • Electrons, at high onorglos 
will bo stooponad by onorgy losses by a factor of E* 1 from thoir source 
spectrum. Electron data above 20 GeV exhibits a spec true as steep as 
E*’ 3 which is consistent with an E 2 * 3 source spectrus. 

Recent studies using HEAO-3 data (Ormes and Prothoroe, 1983) and 
balloon.daua (Simon and Mathis, 1983) now suggest an escape steepening 
of E~* 74 * 1 . These studies imply a source spectrin of E 2 ^which 

in turn suggests an energy steepened electron spectrum of E 3>ofS * 1 , 
which is inconsistent with the electron data. In this paper, the 
authors suggest that a truncated pathlength distribution due to a lack 
of nearby sources could be responsible for the additional steepening of 
the electron spectrum. 

II. Electrons in the Leaky Box Model 

In the leaky box model it is assumed that cosmic rays are Injected 

-r© 

into the galaxy with a power law energy spectrum, N(E) • HI , and 
escape in an energy dependent manner with an average escape lifetime, 

-6 

t_ ■ t E .It is also assumed that these sources are homogeneouwly 
distributed and. that cosmic rays diffuse from these sources throughout 
the galaxy. This suggests an exponential age distribution: 

m ^ ' (z) 

EZ 

Due to cosmic ray electrons* energy losing interactions with the a . 
interstellar electromagnetic fields, electrons have a second cha'rac- 
teristic lifetime. This second lifetime is an energy lose lifetime, 
t ■ l/b£, where b is a constant depending only on the galactic 
electromagnetic fields. 

Using this model, Sllverberg and Ramaty ( 1973 ) and Ramaty (1974) 
have modeled the electron propagation. These authors demonstrate that 
if 6 < 1, at low energy the electron lifetime will be dominated by 

* NAS/NRC Research Associate 
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wupi (t_ « t_), ud at high roarer the electron lifetime will be 
dominated*!? energy Xoaa (t m « t__). Tha equilibria alactron 
apaccrta will ba tha injection spectrum steepened by aaeapa me lav 

-<r fl ♦ «> 

energies (W(I) « * ) and will ba tba lnjaction apaetraa 

. “< r 0 ♦ 1 ) 

steepened by energy loas at high energies (K(£) « E ). 


3. A Truncated Pathlength Distribution 


In raallry, tba distribution of eoanle ray aourcas la not 
uniform. Tba dlatanca no tba naaraat eoanle ray source will datemlne 
sons' nlnlnua llfatlae f«,.r any observed eoanle rays. In tbla work, tha 
authors have used a truneatad age distribution of tha fora: 


P(T)dT 


dT 




whara t_ • t, 4- t_. i g (E) la tba time It takas far a eoanle ray of 
energy fr to diffusa to the solar systen from tba nearest source(s), 
and t la tha average storage Ilf atlas In tba galaxy. If tba energy 
dependence of eoanle ray propagation la assumed to be uni fora throughout 
tba galaxy, than t. and t, are also propolrtlonal to t’ . This age 
distribution is ld£ntleal L to tha age distribution of Cows lk and Wilson's 
(1975) two-sone distribution, cxeapt that t_ new ebaraetarlsaa tba loeal 
storage volume rather than that of the sources. 

Suppose Tg » ft., where f Is some eonstant fraction than: 

T e - ft E* 4 and • (H)r l^. 
so L o 

Considering only synchrotron and Conpton energy losses, the 
equilibrium alactron flux obeys the equation (Silverberg and Eamaty, 
1973): 

If cwVe» *r% ■ «*> 

EX 

-T 0 

with tha source spectrum, Qi’E) a E , the electron flux ran be 
calculated for this truncated pathlength distribution to be: 

-(r +i) -r 

*(E) - W o E 0 ^det ° {exp(-y/(l-f))-exp(-y/f)] 

where 

y - i and c - (l-bE9)~ l 


The variable of Integration, c, is a function dependent on 9, the amount 

of time since the electron left Its source. Sines « E~* and « 

6* S 
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l"*, y m mad y approach** taro aa E Increases, at aoaa wry 

Urge E It la possible to uaa tbs linear approximations: 

s*p(-r/(i-f)) • i - y/(i-f) •** <-r/*> ■ (W '/f) 

Da log those substitutions, ,th* elsetron flux can bo written as: 


H(E) 


." (r o * 




T EH .. /o ,1 . 1 v (If 6 " 1 ) 

tT Ji dcc ( i 


The integral is now independent of E and is a constant, and at eery 
high wnergy we see: 

-<r + 1) -(r ♦ 2 - 8 ) 

N(E) « E %/T n « E 


The injection spectrum has been steepened by a factor of (2 - 6). 

4. Discussion 

The above shows that the lack of nearby sources can steepen the 
injection electron spectrum by more than one power of E. The electron 
spectrua can be broken into three regions. The first is the low-energy 
region where electron storage is dominated by escape, < T„ t . The 
slope of the electron spectrum, like other primary species, is the slope 
of the injection spectrum steepened 6 by escape. The second region is 
the medium energy region where t_ < t_ . < t__ . In this region, the 
electron energy loss lifetime is suf flcientiy short to dominate 
propagation, but not so short as to prevent electrons from propagating 
throughout the storage region. In this second region, the slope of the 
spectrum Is equal to that of the Injection spectrum steepened by one 
power of energy. The third region is very high energy where t < r«. 
This is the energy region In which the energy loss lifetime is shorter 
than the time it takes for cossiic rays to diffuse to Earth from the 
nearest source. This produces a further steepening In the electron 
spectra. 

The salient point of this calculation Is chat the asymptotic slope 
of the electron spectrum can be steeper than that of the injectior, 
spectrum by more than one power of energy. Figure 2 shows the electron 
rpectra above 1 Gc7 as measured by Nishimura et al. (1980) and Prince 
(1979). The steepening of the electron spectra appeals to occur in the 
region of 20-50 GeT. For this steepening to occur due to s truncated 
age distribution, at 50 GeV: 

T s ' ,T n > t em ' ( “>' 1 

The energy loss coefficient, b, has the value 

b - (3.79 x 10" 18 eJ + 1.02 x IO" 16 PpH >(CeV-s)* 1 


<1 


originm- p** ® 


OC3-5 


vhere £Lp 1* the tvirtge trenaveree magnetic field in mlcrogauss end p 
la the photon energy density in eV/em 3 . Using the — 1 *— * 1 
end .7 eV/em 3 for p«„ b - 1.66 * 10“ lc> CeV’ 1 a" 1 . 


Of P, 


values of SwG for 
For this value 

at SO GeV one gets an energy loss lifetime of 3.8 million years. Using 
a value of .2 for f , t„ must be on the order of 20 million years at SO 
GeV. If one assumes tne same storage time rigidity dependence as 
derived by Ormes and Protheroe (1983) for the escape length, then this 
lifetime could be a factor of two longer at 1 GV/c. This can be 
compared with the 10-20 x 10® year lifetime measurements at 1 GV/c based 
on the ^®Be isotope (Wledenbeck and Greiner, 1980). Assuming these 
lifetimes should be the same implies that the root mean square magnetic 
field must be in the range of 8-10uG. 


Although the truncated age distribution provides an attractive 
method for steepening the electron spectrum, it requires that cosmic 
rays are either stored in the galaxy longer than currently believed, or 
the magnetic field traversed by the cosmic rays must be higher than the 
5pG used in the above calculation. 
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Dapt. of Physics, University of Adalalda, Si h Australia 5001 
1. INTRODUCTION 


Considerations of tha Ilf a time, L/K ratio, low anisotropy, and high energy 
alactron spectrum suggest that observed cosmic rays proptjate In a local spaca of 
dimensions lass than 1 kpc. Kafatos at al. (1981) have suggested that tha solar 
systaa may be Inside a super bubble (Bellas, 1979) which confines local eosmlc 
rays. Bruhwailar at al. (1980) and Toaisaks at al. (1981) have shown that tha 
shall-llka structure of auperbubblas may result from tha sequential explosions of 
supernovas In evolving OB associations. In this paper vs explore the consequences 
for cosmic ray phenomena of the solar system being Inside a superbubble which we 
suggest should be Identified with tha expanding ring of BX (Feature A), which was 
first observed by Lindblad (1967). Figure 1 shows tha ring relative to the solar 
system In tha plana of tha galaxy. Figure 2 la a model of the local galactic plana 
magnetic field, which has bean largely excluded from tha suparbubbla by being swept 
up as matter Is "snow-plowed* outward by tha suparbubbla expansion. 

II. LOU ENERGY COMPOSITION AND ENERGY SPECTRA 


The radius of the suparbubbla at age t la proportional to t**(0.5). This 
expansion with time causes tha contained relativistic cosmic rays to lose anargy 
as dE/dc - (-l/2)(E/t). We shall adopt this anargy loss law also for eon-relati- 
vistic particles. This may result In soma underestimation of mafer traversal and 
superbubble age from low energy data. We aaauae that tha Injection rata of nuclei 
of type i of energy E (per nucleon), Q t (E), la constant. Tha region Inside tha 
superbubble la cf low danalty. Cosmic ray particles traverse the Inside region 
relatively freely, interact with the dense walla of the auperbubbla, and accumulate 
grammage In repaated wall encountars throughout their lifetime. On average, over 
repeated wall encounters, the rate of matter traversal by cosmic raya (in g cm" 
par sacond) la taken as Independent of t and paramatarlsad as X(E)/T whera T la tha 
present age of the superbubble. 


The number and energy spectrum of nuclei of type 1 contained vlthln the 
superbubble throughout lta history la given by: 


?N t (E,t) 


Q,(E) ♦- 

T 




Nj(E.t) 


X(E) 

X )1 


( 1 ) 


N (K.t) | 

T x l 




' 4 "5i r V 1- '*' - 


where: x, I* the Interaction langth (g cm - *) of apaclaa t; x,. l » transformation 
length for J ♦ 1; t* (E) la the mean radioactive decay time of species 1. 


Because of the euperbubble 1 s expansion, particles obssrved now with enorgy E 
had energy E’ • E x (t**(-1/2)) when produced at '•poch x I t/T. Solving equation 
(1) for the present spectrum, we obtain: 

l 

N (E, T - l)dR - / P (E\T)dE* ®(E,T) ♦(E.T)dT, 

1 0 


( 2 ) 
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where P.(E' , T)dE' - {Q (E*)T + I N.(E’,r) \ T _l/2 dE. (3) 

1 1 j>i J X J1 

In equation (2), 8(E,t) and 4(E,t) arc ( the probabilltlaa of a nucleus of present 
energy E having survived nuclear Interaction or radioactive decay respectively fron 
production at an earlier epoch, t, with energy E x (t**(-1/2)). 

For a primary species Injected uniformly throughout the history of the super- 
bubble with a power lav energy spectrum with exponent -t» 

P^E'.OdE' - T ^"^QjCEjTdE. (A) 

The quantity analagous to the age distribution of primaries in_eoovent lonal models 
Is t**((y-1)/2) and we will refer to It as such. For X(E) « E , the age distri- 
bution of secondary species Is approximately proportional to r**(l + (y + v ” 
l)/2). These age distributions have been plotted In Figure 3 for y ■ 2.7 where they 
are compared to age distributions In the leaky-box nodel (exponentials) for various 
mean escape times, t e#c * 

Present spectra are similar to production spectra, e.g., primaries at high 
energy (X « x<) have N,(E) - 0.54 Q.(E)T for Y - 2.7 (the 0.54 factor is due to 
expansion). Although trie age distributions of primary and secondary nuclei are 
intrinsically different, secondary and primary energy spectra differ mainly because 
of the assumed energy dependence of X(E). Predicted secondary to primary ratios 
using X * 7 g cm for (E ~ 2 GeV/nuc) and X * 7 (E/2 GeV/nuc)**(-0.33) g cm for 
(E > 2 GeV/nuc) give a good fit to the observations (e.g., de.ta surveyed by 
Protheroe, et el., 1981). Using the observed surviving fraction of l0 Be (Carcia- 
Munoz et al. , 1977) and equations (2)-(3), we obtain T ■ (2.9, +1.3, -0.7) x 10**7 
years, consistent with the astrophyslcal estimate of the age of Feature A (Olano, 
1982). This Is just what we expect based on Figure 3 where we see chat the age 
distribution of young secondaries Is very similar to that In the leaky-box model 
provided T Is about 2.3 times the mean escape time, t e#c , that t uuld be obtained if 
the leaky box model had been used to Interpret the data. 

111. HIGH ENERGY COSMIC RAYS 


The local superbubble model offers a natural explanation for features in the 
high energy cosmic ray anisotropy and spectrum which occur around 10**15 eV (see the 
review of Llnsley, 1981). We interpret these a6 due to failure of the superbubble 
wall to contain cosmic rays of high energy. At energies well below 10**15 eV, the 
wall is assumed to trap completely incident particles. However, at high energies 
the wall becomes more permeable to cosmic rays and there is an energy dependent 
probability G(E) that a particle striking the wall will escape. This is relatec* to 
the high energy anisotropy: 4(E) ■ G(E) Ap(E)/p(E) where p is the cosmic ray 

density and Ap is the difference in density between cosmic rays Inside and outside 
the superbubble. 

At high energies, where we can neglect nuclear interactions, the number of 
cosmic ray nuclei of type i within the superbubble la given by: 

3N.(E,t) . , - . 

—L - Q 1 CE) + + ( Pl (E)V(t) - N 1 (E,t))/tJ;(E,t) (5) 

where E is now taken to be the total energy per nucleus, V is the 1 super bubble's 
volume, P t is the exterior cosmic ray density of species i, and t* le the mean 
leakage time (In or out) of species 1. We shall assume a rigidity degendenee for 
the leakage time of species 1 (atomic number z), t^ (E) - <t>(E/zE Q ) where <t> is 
the mean cosmic ray age (~ 0.4 T) and E is the energy at which the present leakage 
time of protons is equal to <t>. For the exponent of the escape law, we shall use 
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o - 0.6, Che power lew exponent of the observed anisotropy above 10**14, since the 
escape probability la approximately proportional to 6. (We can neglect the energy 
dependence of Ap/p unless the exterior 'spectnn Is steeper than the Interior 
spectrum). Equation (5) nay be solved in two extreme cases for the energy spec trun 
at the present tine (t “ T) to yield: 

Q. (E)T/[1 + (y-1)] (Expansion dominates) 

V E) “ , , (*) 

{p t (E) V Q 1 (E) t*(E)Hl-'exp(-T/t*(E))> (Leakage doninatea) 

For the purpose of illustration, we consider cosoic rays to be composed solely 
of protons and Fe-nuclei. In Figure 3 we show spectra and relative abundances 
calculated using an approximate solution to equation (5), which agrees asympto- 
tically with equation (6), for three representative sets of parameters. In case I, 
the parameters represent minimum assumptions required to obtain agreement with the 
observed total energy spectrum: the spectra produced Inside and outside the 

superbubble are normalized to data at 10**10 eV and 3 x 10**18 eV respectively and 
taken to have the same differential spectral index y • 2.7. E Q has been chosen as 2 
x 10**15 eV and, at a given total energy, a production ratio p:Fe of 3:1 has been 
used. Figure 3(b) shows components of the differential spectrum from interior and 
exterior sources. Case II is a modification of Case I in which a spectral index of 
Y * 3.05 has been used for spectra exterior to the superbubble and a has been set 
equal to 0.8. This causes the spectra from the exterior to dominate at all energies 
above — 10**15 eV and produces a bump-like feature at ~ 10**15eV. Case III is a 
modification of Case I in which Fe is injected with a flatter spcctrim (y ■ 2.6) and 
E q is lowered to 3 x 10**14 eV. As seen in Figure 3(c), this results in an enhanced 
Fe composition consistent with the Maryland result (Goodman et al. , 1979) — (hatched 
area) and the Durham group (Chantler, et al. , 1983) — (vertical bands). 

In the energy range from 3 x 10**14 eV to 10**17 eV, the direction of the 
measured anisotropy (See Linsley, 1981; Watson, 1982) is RA - 12-18 h indicating a 
net local flow from the nearby wall (see Figure 1). Although this might imply a 
higher cosmic ray density at these energies outside the superbubble, this is 
unlikely since the leakage time constants above 10**15 eV become much less than the 
cosmic ray age; a significant net inflow cannot then be sustained. Thus, the 
direction of observed anisotropies from 10**15 eV to 10**17 eV is probably a 
transient of duration much less than <t>. Possibilities Include: a) Cosmic rays 

from exterior sources in the Sco-Cen active region may be "seen" through the partly 
transparent superbubble wall; b) A mini -superbubble (Davelaar, et al. , 1980) may 
surround the Sco-Cen region (Loop I, the North Polar Spur, may be the Intersection 
of the two superbubbles) and contain a higher cosmic ray density which is now 
leaking into our superbubble; c) We are just seeing an anisotropy from the nearest 
(wall) source. Just above 10**17 eV, the abrupt reversal in the observed anisotropy 
direction indicates a return to net outward flow of cosmic rays toward the local 
wall. This persists until about 5 x 10**18 eV, the apparent energy at which the 
anisotropy ceases being dominated by the local wall. 
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ON THE LOCALIZED NATURE OF THE GALACTIC COSMIC RATS 


J, F. ORMES 

NASA/Goddard Spaca Flight Cantar 
Greenbelt, MD 20771 USA 


1. Introduction 

The recent observations of charge composition data from the HEAO-3 
satellite (Engelmann at al., 1981), have bean interpreted (Ormas and 
Protheroe, 1983) to suggest that the escape of cosmic rays from the 
galaxy at rigidities above 10 GV/c falls more rapidly than previously 
thought X a R , (6 • 0.7 ±0.1 compared to earlier values in the range 
0.3 < 6 <*0.5). If one associates this steep rigidity dependence with 
the diffusive escape of particles from the galaxy’, and further connects 
it with the energy dependence of the anisotropy above 10 GV/c, one is 
led to speculate that j^he diffusion coefficient of 10 CV/c cosmic rays 
■ay be as low as 10 26 “~ 7 cm 2 /sec. In this paper we explore the 
consequences of such a low value of the diffusion coefficient and of 
various values of 6. We then address consequences for the propagation 
of cosmic rays, in particular their localisation near their point of 
origin sufficient to explain the roleover the escape of low energies. 

2. The Diffusion Coefficient 

Consider a model in which the diffusion coefficient is isotropic 
throughout all space, and the sources are confined to the galactic disk 
(e.g., where SN are occurring). It can be shown that the spectra 
integrating over the finite disk and for all time give power law spectra 
steepened from the source spectra by the power £ with which the 
diffusion coefficient varies as the rigidity (Streitmatter and 
Balasubrahmanyan, 1983). 

The observed cosmic ray proton spectrum is presumed to be due to 
the injection spectrum, perhaps from shock acceleration, steepended by 
the energy dependence of the diffusive escape of particles from the 
galaxy. Since the proton spectrum has the form of a power lav without 
change of slope from 10 GV/c up to 10 5 CV/c (Gregory et al., 1981), we 
assume that the energy dependence of the escape holds up to at least 
this energy and that it is represented by a diffusion coefficient which 
is also a power law in rigidity. We take for discussion purposes a 
diffusion coefficient as shown in Figure 1, the exponents shown being 
6 • 0.5 and 6 • 0.7 as determined by the energy dependence of the 
escape. The former value comes from the preliminary analysis of balloon 
(e.g., Garcia-Munoi et al., 1981) and HEAO-3 data (Perron et al., 1981) 
whereas £ "0.7 is determined using more complete analysis (Ormes and 
Protheroe, 1983; Simon and Mathis, 1983). The £ * 0.5 curve comes close 
to the lower limit on the scattering length given by setting the 
scattering length equal to the gyro radius, whereas the 6 • 0.7 curve 
stays at least a factor of 100 above that limit, even up to 10 s GV/c. 
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Fig. 1. The assumed variation cf the diffusion 
coefficient in interstellar space. 
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Fig. 2. The saplitude of the anisotropy is well 
fitted by a diffusion and a Compton* 
Getting component. 
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(It it worth noting here thet the escape length may not continue to 
decrease towards higher energy (see Chappell and Webber, 1981), In which 
case the connection between high (>10 GV/c) and low (<10 GV/c) energy 
regimes la not dear)* 

12 20 

The anisotropy amplitude in the range 10 1 to 10 eV Is shown In 
Figure 2. It is well fitted by two components: at energies below 10 14 
eV is a Compton-Getting anisotropy which corresponds to a streaming 
velocity of 20 ka/sec and above is a diffusive anisotropy w|j£c^ g^mes 
from the cosmic rays flowing down the intensity gradient fc e " ’ K ~ d , r w ^ ere 

N is the numer density and r is the spatial coordinate. A. diffusion 
coefficient power law 6 in the range 0.5 to 0.7 fits this data well. 

Using the same diffusion coefficient at 10^ Gv/c and the anisotropy 
of 10” 3 at 10® GV/c, we can compute a half height for the intensity 
gradient of 1.5-15 kpc depending upon the value of 6. These surely 
represent upper limits for the characteristic dimension of the storage 
volume because of the proximity of the solar system to the center of 
symmetry of the galactic disk. 

From the time scale inferred from *®Be observations, we know that 1 
GV/c nuclei have an age of about 10-20x10® yrs. (Wiedenback and Greiner, 
1980 and Garcia-Munoz , Simpson and Wefel, 1981). Using this as a 
characteristic time, we see that with the small diffusion coefficient, 

10 GV/c particles (4 GeV/nucleon) can only diffuse 100-200 pc from their 
sources whereas 10® GV/c particles diffuse 3-10 kpc. It is unlikely 
then that there are more than one or two sources this close to the 
sun. The absence of a break in the high energy electron spectrum has 
been Interpreted as saying that the nearest source must be within a few 
hundred parsec of the sun, but the spectrum of electrons is inconsistent 
with 6 *0.7 and a standard leaky box picture. This problem may be 
resolved if the pathlength distribution is truncated due to the 
separation of the sun from this source. (See accompanying paper, Mauger 
and Ormes, 1983). 

In summary then, we find that at low energies we are 100-300 pc 
from the nearest source, and that at higher energy the diffusive length 
is in the range of 1-10 kpc. These distance scales apply asumlng that 
the diffusive medium has a characteristic size larger than these 
dimensions, and is the same everywhere throughout the volume. If the 
storage volume has smaller dimensions, then physical dimensions will be 
controlling, and the lifetime to reach that boundary will be inversely 
proportional to k . 

3. Conclusions 

The relatively small distance which low energy particles (< 10 
GV/c) can diffuse suggests that the roleover in the matter traversed 
below this energy (Garcia-Munoz et al. , 1979 and Ormes and Protheroe, 
1983) may be due to particles not yet having diffused to us from the 
nearest source. This is consistent with an Increase in the truncation 
of the pathlength distribution at lower energies as reported by the 
Chicago group (Wefel et al- 
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Such * •■all diffffio6°coef9icient 1« consistent with the 
observation by Cana (1977) of 10-20 MHz radio emisslvity on the lint of 
tight to nearby HI1 regions. This data suggests there Is a high 
eaisslvlty region near the solar system. It would also imply the 
localization of particles to create point like gamna-ray sources. A 
prediction of this model would be that gaoma-ray sources would be bigger 
at higher energy. Such an effect might be observable with the extended 
energy range and high sensitivity of the energetic gamma ray experiment 
(Fichtel et al., 1983) on the Gamma Ray Observatory. 

The low energy cosmic rays would then have- been confined very near 
the solar system, so that the observed grammage traversed and the 10 Be 
abundance which lead to the conclusion that these cosmic rays have 
sampled a low density regioi of space p - 0.3 atoms/cm 3 . This low 
density is consistent with the very local region of the galaxy where the 
sun may be inside a hot low density plasma or superbubble blown in the 
cold material of the average interstellar medium by repeated supernova 
explosions. 

This picture leads to the following mix of sources. At low energy 
where the diffusion coefficient is small, each source is relatively 
isolated, and the composition we observe may be unique to that local 
source (and hence its general similarity to the solar system 
material). On the other hand at high energies in which our knowledge 
comes mostly from air shower observations, the diffusion coefficient is 
large, and we are bathed in the cosmic rays produced in a large number 
of sources, hence the low anisotropy values. There are no observations 
on the long term constancy of the cosmic ray intensity which rule out 
fluctuations on a time scale of 10 7 years (Schaeffer 1977). 
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ACCELERATION AND PROPAGATION OP GALACTIC 
COSMIC RATS: IMPLICATIONS OF HEAO-3 DATA 

J.F. Ones 

NASA/ Goddard Space Flight Centra, Greenbelt, Md 20771, U.S.A. 

R .J . Protharoa 

Dapt. of Physics, University of Adelaida, South Australia 5001 

ABSTRACT 

We rc-exaaine the energy dependence of the 
mean escape length of cosmic rays from the galaxy in 
the light of recant measurements of cosmic ray 
abundances from the Danish- French experiment on 
HEAO-3. Implications for cosmic ray acceleration and 
galactic confinement are discussed. 

1 . Introduction 

Highly accurate data on the relative abundances of boron 
through nickel are now available from the Danish- French collaborative 
experiment (Bouffard at al., 1982) on HEAO-3. The boron to carbon ratio 
obtained from this experiment (Eagelmann et al., 1981) is shown in 
Fig. 1 together with low energy observations from the IMP satellite 
(Garcia-Munoz et al., 1979). Balloon observations (original compilation 
by Garcia-Munoz et al., 1981a, with additions) are also shown to 
Indicate the general agreement between these and the more accurate 
satellite data. We have used these results and performed a propagation 
calculation to determine the energy dependence of the mean escape 
length, X (E), from the observed secondary to primary ratios. 
Astrophysical implications of our result will be discussed. Full 
details of our analysis can be found in Ormes and Protheroe (1983) . 

2. Mean escape length 

The methodof calculation is as described by Protheroe et al. 
(1981) . We have used source elemental abundances derived by Perron et 
al. (1982) and assumed source Isotope ratios to be solar (Cameron, 
1980), except for C, 0, Ne, Mg, and Si for which we used those obtained 
by Viedenbeck and Greiner (1981) . For source abundances of the sub-iron 
group we have taken the local galactic abundances (Meyer, 1979). Energy 
dependent total cross sections of Letav et al. (1982) were used. For 
spallat'on cross sections we have used the seml-emplrlcal formulae of 
Silberberg anc Tsao (1977) (see also Tsao and Sllberberg, 1979 and 
references therein) except that those of iron have been normalized to 
measurements of Webber and Brautigas (1982) at 980 MeV/nuc. 

First, we shall restrict our analysis to the high energy data 
(2.6-15 GeV/nuc) to probe the assymptotic behaviour of X , avoiding 
biases introduced at low energies by solar modulation effects, strong 
energy dependence of cross sections, and velocity dependent propagation 
effects. For a source spectrum appropriate to acceleration by strong 
ah.cks, i.e. dJ/dT «p" 2 , we have calculated the secondary to primary 
ratios for two possible forms of the energy variation of X : (a) a power 
lav in rigidity; (b) a power lav in kinetic energy. We shSv in Figure 2 
the results of a comparison of our calculations with three important 
secondary to primary ratios obtained from the HEAO-3 data for the two 
forms of X (the equations used are given in the figure). The results 
sre given in the fora of a x 2 contour plot in the A-d plane. For 


2. 


ORIGINAL page is 
O f POOR QUALITY 


OC 5.2-14 


e <u 
»- 

ae 

« 05 

s 

m 

5 02 
g 

© 01 


h/ 


if GARCIA- MUNOZ ft 01,1979 
• ENGEIMANN eld, 1981 
I BALLOON OBSERVATIONS 


- 

1 \ 


a — uj 


JUUJ 


i,i,A„ixjal 


JJ J-Ui 


0.01 


100 


Figure 1 
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KINETIC ENERGY ( GeV/ NUCLEOK ) 

Boron tc carbon ratio measured by HLAO-3 and IMP satellites 
compared with a survey of balloon observations. 


clarity, only the contour corresponding to the 95! confidence interval 
has been plotted. It is clear that for the B/C, K/0, and Sc-fe/Fe 
ratios, the BEAD data give a best value of 6 somewhat higher then 
previously realised. Vie find 6- • 0.7 ± 0.1 for the rigidity dependent 

fit and 6_ * 0.63 ± 0.1 for the kinetic energy per nucleon dependent 
fit. T 

Previous estimates based on balloon data placed 6_ in the 
range of 0.3-0 .5, however Perron et al. (1981) in th*ir analysis of the 
wvAn data suggested larger values of & would fit better. ‘Tf one were to 
Ignore the quoted errors on the HEAD data and Include data at lower 
energies in the fit, lower values of 6 would be obtained. Also, 

kinematic effects make 6- larger than 6. in this energy range. 

Prom Figure 2, we find that the best value of A obtained from 
the ratio of iron-secondaries to iron is about 10! higher than that 
obtained from the boron to carbon ratio. This is consistent with a 
slightly truncated pathlength distribution, for example as expected in 
the nested leaky box model (Cowslk and tfilson, 1973; and would imply 
about 8! of could be in source regions. Unfortunately, this 
conclusion cannot be reached because of systematic uncertainties it 
spallation cross sections of order 10!. To illustrate this point the 
acceptable range of A-» obtained from the boron to carbon data allowing 
for a 10! uncertainty in the partial erosa sections, has been added to 
Figure 3(a) and Includes tbe region between the two dashed lines. The 
'best' values of A_. obtained from the three secondary to primary ratios 
shown it Figure 3(a) then appear to be entirely consistent. For further 
discussion of truncated pathlength distributions, see Protheroe et al . 
(1961) and Gards-MUnoz et al. (19£la) . 

Vie show in Figure 3 tbe naan escape length we derive frot the 
boron to carbon ratio results from the HEAD- 3 and IMP experiments. For 
this, we have assumed that the solar -modulation cat be described by tbe 
force field approximation with a acceleration parameter, 4, of 600 W>* 
appropriate to near solar maximum conditions (Urch and Gleason, 1973). 
Ve show (dashed lime) 1 • 35 (l/CV/c) * 7 g/cw 2 we obtained for date 

above 2.8 GeV /sue . Clearly, this gives a poor fit to the low energy 

data. Satisfactory a gree m e nt with the data is obtained (solid line) 
with, 
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X, - 3 J ( 1 ♦ P (^)*« •’«/«? (1) • 

The turnover in A # at low energies In considerably more abrupt than ean 
ba accounted for bp the Dynamical Halo model (see e.g. Jones, 1979). 

3. Pis cuss Ion 

At high energies, the variation of A € vith energy msy be as 
steep as *T® * 7 * Since the observed spectra of primary nuclei (e.g. lyan 
at al. 1972) are proportional to IT 2 * 7 , this would Imply that primary 
cosmic rays may be produced with a p" 2 spectrum as expected for 
acceleration by s first order Fermi mechanise In strong shocks 

(compression ratio 4, or certainly greater than 3.5; see e.g. Axford, 
1981). 

As well as implications for the aeceleradon of cosmic rays, 
the steep rigidity depender.ee of A e has important consequences regarding 
the distance cosmic rays propagate (the cosmic ray age must be greater 
than the speed of light crossing time of the storage volume). Since 
there is no observed structure (e.g., change of slope) in the proton 
spectrum up to 10 s GV/c (Gregory et al., 1981, Tasaka et al., 1982), we 
conclude that the rapid decrease of A continues up to this rigidity. 
Vie shall assume that the lifetime of cosmic rays is proportional to A . 
Extrapolating from «*10 7 y at 1 GV/c (Wiedenbeck and Greiner, 1980; 

Garcia-Munoz et al., 1981b) using equation (1), the age at 10 s GV/c is 

4000 y sod Che size of the "storage region" must be less than 1 kpc . 
Particles which propagate diffusively must strongly satisfy this 
inequality. 



Figure 2 Comparison of boron to carbon, sub-iron to iron, and 
nitrogen to oxygen ratios observed from 2 .8-15 GeV/nuc by 
HEAD-3 with results of the present propagation 
calculations . The figure shows the goodness of fit 
when A. end depart a) or A. and 6. (part b) are varied. 
The cl Seed curves ere contours of constant x 2 corresponding 
to 952 probability. 
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fart of this work was performed as a KEAO-3 gusst 
Investigation. J.7. Oraas would Ilka to thank Drs. L. koeh-Mlraaond and 
I. Lund for thalr hospitality and discussions during visits to Saclay 
and DfXZ in 1981. We art grataful to Drs* 8 *A. Stephens, 
?4C. Belasobrahnanyan and M. Slnon for useful discussions* 
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Figure 3 Mean escape length derived by us fron HEAO-3 and Itt 
aeasuraaents of the boron/carbon ratio. Dashed curve shows 
assyaptotlc behaviour; solid line la the fit given in 
aquation (1). 
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HOW DO COSMIC RAYS CHANCE THEIR ENERCY IN THE SOLAR WIND? 

Frank C. Jonas 

Laboratory for High Energy Astrophysics 
NASA/Goddard Space Flight Center 
Greenbelt, Maryland 20771 
USA 

ABSTRACT 

We have derived the diffusion-convection (modulation) 
equation directly fros the Boltzman equation employing a 
minimum number of assumptions about the scattering 
process. We assume that each acatterer has s reference 
frame in which 1) the scattered particles undergo no change 
in energy, and 2) isotropy Is an equilibrium state. If the 
background plasma contains a magnetic field and the flow 
speeds of the plasma and scattering centers are different 
additional terms arise that modify the usual equations. 
Furthermore if the scatterers have Individual motions 
relative to their average flow the second order Fermi 
acceleration term appears. 

1. Introduction 


The dif fusion-convection equation was first applied to cosmic-ray 
modulation by parker (1965,1966). Gleeson and Axford (1967) were the 
first to derive It from the more general Boltzmann equation: 



♦ 

v 


$ f + ? 


$ f 
P 


-1 
3t c 


( 1 ) 


where the left hand side of the equation represtr.-ts the streaming of the 
particles through the large scale, average force fields represented 
by F and the rlghc hand side represents the change in the distribution 
function f due to stochastic forces or scattering. In their paper 
Gleeson and Axford did not include any large scale magnetic or electric 
fields and treated the scattering as hard sphere scattering. If one 
considers a magnetic field embedded in a plasma flowing with a velocity 

v one must Include both the magnetic and induced electric fields In the 
P 

force term, 

f - - * — - m V p x w - (?-my tf p ) * (^) (2) 

where p is the particle momentum, m Is the particle mass, Y is the 
particles Lorentz factor, V p is the plasms flow velocity 

and a! » ei/mc, the particles local gyrofrequency . 
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2. The Collision Tern 

We may calculate the colllalon tent In the aanner of Gleeaon and 
Axford by performing the calculation In the rest frame of an Individual 
scattering center. In this frane ve sake two important assumptions 
namely: the scattered particles do not change their energy in a single 
scattering and the scattering process is such that it will preserve 
Isotropy. We further assume that the momentum transfer, averaged over 
the cross section, may be written as a tensor operating on the original 
momentum. 

It should be noted that ve assume nothing about the strength of the 
scattering; it could be strong as in hard sphere scattering or weak as 
in scattering by weak Alfven waves. 

3. The Transport Equation 

The transport equation is obtained by the usual weakly anisotropic 
expansion, it is assumed that the distribution function may be 
described by its first two moments: 

f(p, x) - n(p,x) + l(p,x)*P (3) 

where n(p,x) ■ Jdpf(p,x) 

and $(p,x) - Jdp p f(p,x). (A) 

If the scattering centers have a mean velocity $ and have a 
dispersion about this mean <6$ s > we obtain after a grSat deal of 

coordinate transformations, integrations and algebra the following 
transport equation: 

- $ • it • tn - I ($ • t p ) p fS ♦ * p . in 

- * 3 pH (5 + v) v • a? p 4».v(iW> (|^n) 

1 > •» «• m —1 r w w m — 1 « ♦ 

+ j AV • v • (ft + v) • [4 - v • (ft + v) I • v n 

+ • v • (ft + v) 1 • ft ^ ^ (mY|>0) 

- d . 5 i + J) . (5 + 5)' 1 . ; . (5 + Jf l . n • |2 

p 8 op 

5 p > If (* y P 3 fj’ < 5 > 


In the above equation K Is the diffusion tensor given by 
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where v in the scattering tensor who** element* arc the various 
scattering frequencies and !) (* the gyro frequency tensor defined 
hy i) • v <o v V /y. It should be noted th*t Inclusion of Q t n the 
d 1 1 fusion tensor produces gradient and curvature drift* of the 
particles. AV I* the different** between the plasma *nd acattering 

root or velocities, AV V - V and the last term on the right hand aid* 

P « 

t* tho Form! aooolorat Ion torm and I* *ero If the scattering centers *11 
mo vo with a single voloolty, l.o., AV^ * 0. 

4. Conclusions 

two thing* should bo noted about equation (A): flrat, if the^ 
aoattorlng frequency t* small compared to the gyro-frequency (\><<fl), a* 
Is usually tho case, all of tho terms on the right hand side (RHS) are 
small and mav bo neglected (with the possible exception of the Fermt 
terra); second, tho velocity that appears In the left hand side (UIS) 

(tho usual modulation equation) Is the plasma velocity not the 
scattering center velocity. If we consider the reverse limit, (K<v the 
first throe term* on the RHS are of order uni tv and have the effect of 

replacing V with V on tho l.HS. Tito remaining terms on the RHS are 

P ** 

small in this limit too and arc only of Interest In the Intermediate 
region w • n - 0 and the scattering center motion Is the approriate one 
to use. 


In tho directions perpendicular to ^ , however, the modulation 
equat ton should use the plasma velocity and. In most cases. It Is not 
the mot Ion ot the scattering centers that sweeps the cosmic rays out of 
the solar cavltv but the F \ ft drift Instead. 
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SOLAR ^HE-RICH EVENTS OBSERVED ON I SEE- 3 

D. V. Eenet, and T. T. von Rosenvlnge 
NASA/Goddard Space Flight Center, Greenbelt, MD 20771, O.S.A. 

ABSTRACT 


A scan has been nade of 1SEE-3 data for all 
events with ^He/^He > 0.20 at > 1.3 MeV/nucl. The 67 
events found show soise evidence of 27-day recurrence. 
Larger events show both velocity dispersion and 
magnetic field-aligned arrival from the solar 
direction. At least one third of the events are 
preceded by increases in — 300 keV electrons, although 
several larger events show no electron increases. "^He 
spike events algo exist suggesting nearly scatter-free 
propagation of ’'He from well-connected events. 


1. Introduction 


Despite a decade of observation (see Ramaty et al. , 1980), our 
understanding of ^He-rich solar particle events and of the conditions 
under which they occur remains qualitative at best. In this paper we 
describe some results of a survey of approximately 4 years of ISEE-3 
data that seeks to characterize these events more completely. 

Earlier papers (von Rosenvinge and Reames, 1979; Reames and von 
Rosenvinge, 1981) described the ISEE-3 very low-energy telescope (VLET) 
and showed examples of the element and isotope resolution it provides. 

2. Data Analyses and Results 

The survey was performed using the 1.3-1. 6 MeV/nucl and 2. 2-3.1 
MeV/nucl ^He and *He fluxes in 6-hour averaged intervals from August 15, 
1978 through June 10, 1982. A ■'He-rich interval was defined when either 
or both of the energy intervals met the following two criteria: 

1) The error in the ^He flux from all sources was less than 50 
percent 

2) The %e/^He ratio was > 0.20 (when the Sie flux was zero the' 
ratio was defined as the number of J He particles in the interval). 

A candidate He event consisted of at least two successive Re-rich 
intervals. These events were then observed at a higher time resolution 
to identify multiple events, to account for data-gap effects and to 
better define the onset times. 

Of 5,628 intervals with data, 300 were %e-rich and 67 event 
periods were identified. This amounts to a 5.3 percent probability of a 


2 


ORIGINAL PAGE IS 

OF POOR QUALITY 


SP-2 


%t«rieh event In progress or an average event rate of one every 3 weeks 
during this period of the solar cycle* 

The time distribution of the events Is shown on a 27-day Bartels- 
rotatlon plot In Figure 1* Incur rences seen In the figure exceed those 
expected f roe a randoe distribution and suggest the recurrence of 
events from the sane solar active legion* It has not yet been possible 
to associate eany of * v e snail 3 He events with solar flares to establish 
their origin* 




0.1 I 10 


Fig. 1 27-day 

recurrence 
plot of 3 He-rich 
event times* 


Fig. 2 Scatter plot of 
event -averaged 
3 Re spectral 
Index versus 
3 He/ 4 He 
ratio. 
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OEC 14-15,1979 


Tine histories 
of 3 He sod 4 He 
and Be sectored 
rate plots* The 
direction of the 
Sun is at the 
top and A and B 
are the aniso 
tropy and mag** 
netic field 
directions. 
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Another new feature found in the purvey dete and ahown in Figure 2 
la the tendency of events with large ^He/*He ratios to have steeper 
spectra* 

3. Tine Histories 

A different aspect of the ^He-event characterisation nay be seen 
fron the tine histories as shown in Figure 3 for the largest event in 
the survey. This event typifies two frequently, observed features of 3 He 
events, namely: 1) velocity dispersion (earlier arrival of higher 

velocity particles) and 2) large and persistent f leld-allgned aniso- 
tropies. A third feature of the events Is an absence of correlation 
with interplanetary shocks, sector boundaries or ^He abundance in the 
solar wind. These features establish that the events are Indeed solar 
in origin. 

Further evidence of velocity dispersion is shown in about one-third 
of the events that have clearly associated increases of > 300 keV 
electrons preceding the event by 1 to 12 hours. Electron increases are 
not observed in all events however. 

Figure 4 shows an extreme example of fast spike Increases in ^He 
that are seen In some events. The Increase at 0 hours on November 10 
shows strongly field-aligned flow of 3 He. R. Zwickl (1983) has informed 
us that a "Strahl" Involving field-aligned electron flows described by 
Rosenbauer (1976) is in progress during this period. Efficient 
"scatter-free" propagation of 3 He from an impulsive solar event along 
well-connected field lines could explain the fast profile observed. 
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SOLAR FLARE NEUTRONS AND GAMMA RAY LINES 


R. E. Llngenfelter 

University of California, San Diego, CA, 92093, U.S.A. 

R. Ramaty and R. J. Murphy* 

NASA/Goddard Space Flight Center, Greenbelt, MD 20771, U.S.A. 

B. Kozlovsky 

Tel Aviv University, Israel 
ABSTRACT 

We have derived the energy spectrum of accelerated pro- 
tons and nuclei at the site of the June 21, 1980 limb flare by 
a new technique, using observations of the time -dependent flux 
of high-energy neutrons at the Earth. We find that this 
energy spectrum Is very similar to the energy spectra of 7 
disk flares for which the accelerated particle spectra have 
been previously derived using observations of 4-7 MeV to 2.223 
MeV fluence ratios. 

1. INTRODUCTION 


High energy neutrons, produced by nuclear reactions of flare- 
accelerated protons and nuclei, were recently discovered from solar 
flares with detectors on the SMM (Chupp et al. , 1982). The observed time 
dependence of the neutron flux at Earth gives a time -of -flight measure- 
ment of the neutron energy spectrum at the Sun, provided that the dura- 
tion of the flare is much shorter than the typical Sun-Earth neutron 
transit time (~ 10 minutes). Calculations of neutron energy spectra 
produced In solar flares and predictions of the expected time depen- 
dences of the neutron flux at Earth were made previously (Llngenfelter 
et al., 1965; Llngenfelter and Ramaty, 1967). The time profile of the 
neutron flux observed (Chupp et al. , 1982.; E. L. Chupp, private communi- 
cation 1983) from the June 21, 1980 flare is consistent with these pre- 
dictions. Comparison of the observed time dependence with such calcula- 
tions allows us to determine the energy spectrum and total number of the 
accelerated particles at the site of the flare. 

Neutron production in solar flares Is accompanied by the production 
of nuclear gamma radiation (Llngenfelter and Ramaty, 1967; Ramaty, 
Kozlovsky and Llngenfelter, 1975). Comparisons (Raaaty, 1983) of the 
relative intensities of gamma-ray lines from flares (e.g., Chupp, 1982) 
with calculated line ratios provide independent information on the 
number and spectrum of flare-accelerated particles. 

Comparisons of gamma-ray observations from flares and charged- 
particle measurements In the interplanetary medium have shown (e.g., von 
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Rosenvlnge, Ramaty and Reames, 1981) that gamma-ray production In aolar 
flares takes place predominantly during the slowing -down of the acceler- 
ated particles In the solar atmosphere, rather than during their accel- 
eration or escape from the Sun. This Is reasonable since particle 
acceleration generally requires a law ambient density so that the accel- 
eration rate exceeds the energy-loss rate, whereas effective nuclear 
Interactions require an ambient density that Is high enough to stop the 
particles. Thus, we expect the region of most efficient acceleration to 
be above the Interaction region where the observed neutrons and gamma 
rays are produced. This model of nuclear Interactions Is commonly 
referred to as the thick-target model. 

We have recently (Ramaty et al., 1983) carried out new thick-target 
calculations of neutron production In solar flares. Here we briefly 
describe these calculations; we compare them to the high-energy neutron 
observations of the June 21, 1980 limb flare; and we derive the energy 
spectrum and number of the protons and nuclei accelerated In this flare. 
Then by deriving the fluences of the accompanying nuclear gamma radia- 
tion, we confirm that, for limb flares, the 2.223 MeV line from neutron 
capture on is strongly attenuated by Compton scattering In the photo- 
sphere (Wang and Ramaty, 1974; Chupp, 1982). 

The pectrum of the accelerated protons and nuclei has also been 
derived fvom observations of the 4-7MeV-to-2.223MeV flue nee ratio for 
disk flares. We find that the proton spectra obtained previously 
(Ramaty 1983) for 7 disk flares using this method are very similar to 
the proton spectrum derived using high energy neutron observations for 
the June 21, 1980 limb flare, showing that flares which produce detec- 
table high energy neutron fluxes do not appear to be significantly 
different from other flares observed on the disk. 

II. PRODUCTION OP HIGH ENERGY NEUTRONS 


The calculations presented here are performed using the thick- 
target model. We take the energy spectrum of the protons and nuclei 
that emerge from the acceleration region and are Incident on the Inter- 
action region to be the Bessel function appropriate for stochastic Fermi 
acceleration, N(E) * K2l2(3p/(mcaT) ) l ' 2 ] (Ramaty 1979). The product oT 
characterizes the particle spectrum, such that a larger value of aT 
corresponds to a harder spectrum. The above form for N(E) is valid for 
nonrelatlvlstlc energies and an energy-independent aT. We assume, in 
addition, that the composition of the accelerated particles is the same 
as that of the solar atmosphere and that the protons and nuclei have the 
same energy spectrum. These particles are then allowed to slow down In 
the interaction region due to ionization losses in a neutral medium. For 
simplicity, we also assume that within the interaction region the angu- 
lar distribution of the charged particles is isotropic. An anisotropic 
distribution would have observable effects. For example, if the protons 
and nuclei were preferentially directed downward toward the photosphere, 
the neutron flux from a limb flare would be larger than from a disk 
flare. There is, however, no data as yet to suggest such an anisotropy. 
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Figure i. Observed and calcu- 
lated neutron fluxes at Earth 
as functions of the difference 
between the Sun-Earth transit 
times of neutrons and photons; 
we take d/c-500 sec. 


The calculated time -dependent 
neutron fluxes at the Earth are shown 
in Figure 1 for various values of aT, 
assuming Instantaneous production and 
free, unattenuated escape from the 
Sun. The time At is that In excess of 
the light travel time from the Sun, 
where t Is measured from the neutron 
production time. Also shorn In this 

figure Is the observed neutron flux 
from the June 21, 1980 flare (Chupp et 
al. 1982; E. L. Chupp, private commu- 
nication 1983). We assumed that the 
neutron production time was at 1:18:55 
UT - d/c, i.e. at the midpoint between 
the two impulsive photon emission 
peaks whose centers, shown by arrows, 
were observed at 1:18:40 and 1:19:10 
UT (e.g. Chupp et al. 1982). The 
total duration of the impulsive phase 
of the June 21, 1980 flare was ~ 60 
sec. As can be seen, the shape of the 
observed time dependence fits very 
well that calculated for aT-0.02. The 
absolute normalization of the data to 
the calculations implies a total pro- 
duction of 2.8 x 10 30 neutrons and a 
total number of accelerated protons 
with energies greater than 30 MeV, 

N (>30MeV) * 1 . 2xl0 3 3 . 

P 


The assumption of free neutron escape limits the maximum ambient 
density in the interaction region. If the dependence on height of the 
ambient density in the interaction region can be approximated by an 
exponential with scale height h Q , n(h) « exp(-h/h Q ), then for a limb 
flare, the observable neutron flux originating at a height h Is attenu- 
ated by ~ exp( -( (n/2)Rh ) * n(h)], where R is the solar radius and 

a is the neutron-proton elastic scattering cross section. Since c n p 
increases with decreasing neutron energy } the maximum value of n(h) con- 
sistent with free neutron escape is determined by the lowest observed 
neutron energy A ~50MeV (Chupp et al. 1982). For h Q = 10 7 cm, we find that 
n < 5xl0 15 cm 3 , essentially the top of the photosphere. Since the 
vertical column depth at this density is only ~ 0.2 gm/cm 2 (Vemazza, 
Avrett and Looser 1981), while the stopping ranges of the protons that 
produce the neutrons are >5 gm/cm , the protons must be stopped at 
column depths significantly less than their ranges, perhaps by magnetic 
mirroring (Zweibel and Haber 1983) or by scattering from magnetic 
inhomogene it les . 


Lastly, we consider the 2.223 MeV gamma-ray line resulting from the 
production of neutrons in the June 21, 1980 flare. Calculations (Wang 
and Raraaty 1974) of neutron slowing down and capture on predict that 
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there will be a strong limb darkening of the resultant 2.223 MaV line 
owing to Compton scattering In the photosphere. Taking the estimated 
(Ramaty 1983) value of 0.23 for the number of observable 2.223 KeV 
photons per neutron from a disk-centered flare, the total production of 
2.8xl0 30 neutrons Implies a 2.223 MeV line fluence of ~230 photons/cm 2 
at the Earth from such a flare. The fact that the observed (Chupp 1982; 
D. Forest, private communication 1982) 2.223 MeV line fluence for the 
June 21, 1980 limb flare was only ~6 photons/cm 2 , provides clear evi- 
dence for the limb-darkening of this line. 
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THE GSFC ADVANCED COMPTON TELESCOPE (ACT) 

R. Hartman, C. Fichtel, 

D. Knlffen, G. Stacy' 1 ', J. Trombka 
NASA/Goddard Space Flight Center 
Greenbelt, MD 20771 


ABSTRACT 

A new telescope is being developed at GSFC for the 
study of point sources of gaaaa rays In the energy 
range 1-30 MeV. Using the detection principle of a 
Compton scatter in a 2.5 cm thick Nal(Tl) detector 
followed by absorption in a 15 ca thick Nal(Tt) 
detector, the telescope uses a rocking colllaator for 
fleld-of-vlew reduction and background subtraction. 

Background reduction techniques Include lead-plastic 
scintillator shielding, pulse shape discrimination and 
Anger camera operation in both Nal detectors, as well 
as a time-of-f light measurement between them. The 
Instrument configuration and status will be described. 

1 . Introduction 

There is particular interest in the study of "medium energy” cosmic 
Y rays for a variety of reasons. This energy regime (from a few hundred 
keV to approximately 30 MeV) constitutes a transition region from rela- 
tively low energy thermal processes to highly energetic non-thermal 
astrophyslcal phenomena (e.g., nuclear radiation, s° decays. Inverse 
Compton). For a variety of both galactic and extragalactlc objects 
(supernovae, neutron stars, active galaxies), various models predict 
peak power of emitted radiation In the medium energy Y~ray region. In 
addition, observations of the diffuse Y~ray background emission appear 
to Indicate that there may be Inflections in the energy spectrum of that 
emission In the few MeV region. Such a break In the energy spectrum of 
the y-ray background could contain important cosmological lnforsutlon on 
the history and evolution of the Universe. 

To date, there have been few confirmed detections In the medium 
energy y-ray region. Most experiments have yielded only upper limits. 
The major problem is that in this energy region one encounters low 
fluxes and relatively high levels of Instrumental background; in addi- 
tion, separation of point sources from the diffuse background Is diffi- 
cult because of the limitations of the Y~ray detection process. Below a 
few KeV in energy, there Is a high background due to continuum emission 
and nuclear decay line, arising primarily from cosmic ray Interactions 
In and around the detector. Above a few MeV the background Is somewhat 
less, but the photon fluxes are much weaker, still presenting a sizeable 
slgnal-to-nolse problem. 
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Th*r*for«, th* key to successful oHarvatlom ol aedlua energy y 
ray* 1* th* suppression of Inatruaenta* background. A new Coapton 
t«l*acop* la being developed at GSPC for point aourc* observations 
b«tv**n 1 and 30 M»V. In thla instrument, a coablnatlon of techniques 
la ua*d to r«duc* lnatruaental background. Furtheraore, it la d*algn*d 
to allow accurate aubtractlon of th* real dual inatruaant background aa 
wall aa th* dlffua* ataoapherlc and coaalc flux* a. 

2 . Imstruaent Description 


Th* telescope, which will b* an oriented balloon payload, la shown 
sch*aatically In Fig. 1, Th* central detector* ar* both Nal scintilla- 
tors, each with 1600 cn z ar*a; th* upper detector, in which th* Initial 
Coapton acatter of a y-ray la to occur, la about 2.5 ca thick. The 


lower detector. In which th* acattered gi 
ca in thlckneaa. Each of th* two Nal 
detectors Is viewed by a 5 * 5 array of 
high speed 2" photoaultlpllers . 
Signals froa these tubes are separately 
digitised to provide lnforaatlon on the 
Interaction positions within th* two 
detectors, and therefore the trajectory 
of the scattered y ray. 

The upper Nal detector Is sand** 
wlched between two 6 aa sheets of 
plastic scintillator. A phoswlch 
technique Is used to separate the Nal 
and plastic scintillator signals; th* 
plastic signal is then used as a veto 
to prevent triggering on events In 
which a charged particle enters or 
leaves the upper Nal detector. 

The two central detectors are ' 
surrounded by a shield of lead, about 
5.5 rt thick, which Is encased In 
plastic scintillator anticoincidence 
counters. These are viewed by an array 
of 12.6 c.* photoaultlpllers, 12 at the 
top and 4 at the bottoa. 

Two techniques are esployed to 
ellalnate background froa fast neu- 
trons; In each of the central Nal 
detectors, the pulse shape is deter- 
alned by direct digitisation at a rate 
of 50-70 MHs . This lnforaatlon la 
teleaetered to the ground for later 
analysis. The heavily Ionising nuclei 
resulting froa fast neutron Interac- 
tions produce pulse shapes In Nal which 
are significantly different froa those 


ray la absorbed, la about 15 



Fig. 1 - Instruaent Schematic 

N, Nal scintillators; P, phos- 
wlch plastic anticoincidence 
scintillators; A, anticoinci- 
dence plastic scintillators; 
C, lead alot colliaator; L, 
lead shielding 
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generated by electrons. Th* second t*chnlqu* for foot neutron rejection 
i* a tlme-of -flight e***ureo*nt between the two Wal detectors* With the 
40 ce spacing between these detectors end the relatively slow risetime 
of the Nal pulse, it probably will not be possible to dlrtlngulsh 
between upward and downward aoving y-rays; but it will be possible to 
separate out and reject those events due to the fast neutrons, which 
have typical velocities of c/7 or less* 

The direction determination is made not by the Anger camera Compton 
event reconstruction, but by means of two crossed slat leal collimators, 
giving a f lald-of-view of about 2*. This permits us to use the Anger 
reconstruction as a cross-check on the collimator direction* By 
requiring that the Compton reconstructed direction cone Include the 
collimator direction, it should be possible to reduce by a factor of at 
least ten the backgrounds due to 8-y decays in the Wat or y rays 
originating within the rest of the instrument* 

Despite the background elimination precautions mentioned above, 
some instrumental background will still remain* Furthermore, atmo- 
spheric and cosmic diffuse fluxes must be accounted for in the search 
for localized y-ray sources. This background subtraction will be 
accomplished by comparing the rate of acceptable event while pointing 
on-source to that while pointing off-source. Rather than reorient the 
entire telescope for this purpose, one of the two slat collimators will 
be rocked by a few degrees every few minutes* The on-source and off- 
source positions of the collimator are symmetric with respect to the 
telescope axis, eliminating any bias due to a changing configuration* 

The telescope orientation will be determined and controlled by an 
on-board closed-loop microprocessor-baaed system. This system uses as 
its detector a CID solid state video camera' 2 ' fitted with a 58 am f/1.2 
lens which has extended IR coating. This system has a f ield-of-view of 
8.7* x 11 .2", with angular resolution of 0.045* x 0.036* (244 x 248 
pixels). Preliminary tests indicated that this system can detector G- 
type stars of approximately magnitude 6 with a 5s exposure. 

3 . Conclusions 

Baaed on estimated background as well as Increased effective area, 
it anticipated that, with the instrument described, point source detec- 
tion sensitivity should be Improved by nearly an order of magnitude over 
earlier balloon Instruments. 
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ECRET: THE HIGH ENERGY GAMMA RAY TELESCOPE 

FOR NASA's GAMMA RAY OBSERVATORY 

C.E. Flchttl, D.L. Bertsch, R.C. Hartman, D.A. Knlffen, D*J« Thompson 
NASA/Goddard Spac* Flight Center 
Creenbelt, MD 20771 U.S*A* 

R. Hofstadter, E.B. Hughs*, L.E. Campbell -FI naan 
Hansen Laboratory, Stanford University 
Stanford, CA 94305 U.S.A. 

K. Plnkau, H. Mayer-Hasselwander, G* Kanbach, H. Rothensel, M. Sommer 
Max-Planck-Institflt fGr Extraterrastrlsche Physlk 
8046 Garchlng b. Mtlnchen, Federal Republic of Germany 

A.J. Favale, E.J. Schneid 
Grumman Aerospace Corporation 
Bethpage, NY 11714 U.S.A. 

ABSTRACT 

The EGRET high energy y-ray telescope under development for NASA's Gamma 
Ray Observatory vlll have an energy range of approximately 12 to 30,000 
MeV, energy resolution of about 152 FWHM over most of that range, an 
effective area of about 2000 cm 2 at high energies, and single photon 
angular accuracy of ~ 2* at 100 MeV, < 0.1* above 5 GeV. This instrument 
can locate strong sources to an accuracy of about 5 arc min. The iastru- 
ment utilises a set of digital spark chambers Interleaved with tantalum 
foils for detection and Identification of y-ray events, and a large 
Nal(Ti) scintillator for energy determination. The system is triggered 
by a coincidence matrix using two arrays of plastic scintillation 
counters and a large plastic scintillator anticoincidence dome that 
rejects Incident charged particles. 

1 . Introduction 

The Gamma Ray Observatory (GR0), a Shuttle-launched free-flying 
satellite, is scheduled for launch in the late 1980' s. This will be the 
first mission dedicated to a comprehensive study of the y-ray portion of 
the electromagnetic spectrum encompassing over five orders of magnitude 
In energy. Four Instruments will be included to cover the broad range of 
energies and to achieve the scientific objectives. 

The purpose of this paper Is to describe the high energy Instrument 
to be Included on GRC. This Instrument, a spark chamber telescope with a 
Nal(Tt) total energy counter, will cover the energy range from about. 12 
to 30,000 MeV, with much better sensitivity and spectral and spatial 
resolution than any previous experiment In this energy range. 

2 . Scientific Objectives 

The specific objectives of this experiment are: 

e To search for localised y-ray sources in the energy range from 12 MeV 
to 30 GeV and to measure their Intensity, energy spectrum, position, 
and possible time variations, including periodic emission 
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• To improv* ‘.ha knowledge of the locations of known high energy y-ray 
sources not .identified with objects seen at other frequencies 

• To examln supernova remnants and search for evidence of cosmic-ray 
particle av.-.-. -ration revealed by high energy y r»y« 

• To obtain a detailed picture of the diffuse high energy y-ray emis- 
sion from the Galaxy which in turn will provide information of funda- 
mental importance in the study of the dynamics of the Galaxy. This 
should provide a high contrast picture of galactic structure, inclu- 
ding spiral arm segments, clouds, the galactic center, and the 
extended disk (or "halo”) 

• To determine rhe relative importance of cosmic ray electro. a and cos- 
mic ray nuclernn throughout the Galaxy by studying the y-ray spectrum 
above 12 MeV and, by combining the results with the continuum radio 
data, to obtain a better estimate of the Ralactic magnetic fields 

• To detect and examine the high energy y-ray emission from other 
galaxies, including both normal galaxies and galaxies which are 
exceptionally luminous at other wavelengths 

• To study the diffuse celestial radiation in the energy range above 12 
MeV including its energy spectrum and the degree of uniformity, both 
on a small scale and a broad scale, measurements of particular impor- 
tance in relation to several cosmological models 

• To search for high energy y-ray bursts, to study the spectra of low 
energy y-ray bursts with high statistical accuracy, and to observe 
y-ray lines in solar flares 

3. Instrument Description 

The telescope is shown schematically in Figure 1. Ay ray entering 
the telescope within the acceptance angle has a reasonable probability of 
converting into an electron-positron pair in one of the thin Ta foils 
between the spark chambers in the upper portion of the telescope. If at 
least one particle of the pair is detected by the directional time-of- 
flight coincidence system as a downward moving particle and there is no 
signal in the large anticoincidence scintillator surrounding the upper 
portion of the telescope, the spark chamber track imaging system is trig- 
gered, providing a digital picture of the y-ray event, and the analysis 
of the energy signal from the Nal(Tl) crystal is initiated. Incident 
charged particles are rejected by the anticoincidence dome and low energy 
backward -moving charged particles that do not reach the anticoincidence 
dome are rejected by the time-of-f light measurement. Events other than 
the desired y rays, such as those y rays interacting in the thin pressure 
vessel inside the anticoincidence scintillator, are rejected in the sub- 
sequent data analysis by inspection of the spark chamber images. 

The directional scintillator telescope consists of two levels of a 
four by four tile array with selected elements in each array in a time- 
of-f light coincidence. The upper spark chamber assembly consists of 28 
spark chamber modules interleaved with 27 0.02 radiation length plates in 
which the y ray may convert into an electron pair and in which the inter- 
action point and initial direction of the electrons may be determined. 
The lower spark chamber assembly, between the two time-of-f light scintil- 
lator planea, allows the particle trajectories to be followed, provides 
further Information on the division of energy between the electron and 
positron, and makes possible the observation of the two particles for 
very high energy y rays. 
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The energy of the 
y ray will be deter- 
mined primarily from 
measurements In an 
eight radiation 
length, 76 cm x 76 cm 
square Nal(Tl) scin- 
tillation crystal 
below the lower plas- 
tic scintillator 
plane. Below 200 
MeV, energy informa- 
tion is also obtained 
from measurements of 
multiple scattering 
in the tantalum foil. 
The energy resolution 
of the proposed expe- 
riment is better than 
15Z (FWHM) above 30 
MeV. 


Gamma-ray bursts 
will be recorded in 
both the anticoinci- 
dence dome and the 
large Nal(Ti) cry- 
stal, where an energy 
spectrum will also be 
measured in the 0.3 
to 24 MeV interval. 



The principal Figure 1: High Energy Gamma-ray Telescope 

capabilities of the (EGRET) 

instrument are sum- 
marized in Table 1. For additional details see Hughes et al. (1979). 

4. Capabilities of the Experiment 

The high energy instrument on GRO is characterized by the very low 
background available with the detection of the pair production interac- 
tion. The instrument also has good spectral and spatial resolution and 
sensitivity which should combine to give this telescope the abil-^y to 
locate sources with one-tenth the intensity of the Crab pulsar to an 
accuracy of four to ten arc minutes, depending on the general background 
level. This gives a factor of well over two orders of magnitude improve- 
ment over existing error boxes. The .1 millisecond timing will allow 
searches for pulsars with periods down to the the 1.5 msec range, pro- 
vided the period is known. 

For the diffuse galactic plane emission, the spectra will be mea- 
sured with high accuracy, and spatial variations in the spectrum 
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Table 1 

HIGH ENERGY GAMMA-RAY TELESCOPE 
CAPABILITIES AND REQUIREMENTS 

I. Detectors 

a. Type 

b. Energy Range 

c. Energy Resolution 

d. Total Detector Area 

e. Area Efficiency Factor 

f. Threshold Point Source 

Sensitivity (10 6 sec) 

g. Angular Resolution 


h. Source Resolution 
(strong source) 

1. Field-of-View 

j . Timing Accuracy 

k. Design Lifetime 

2. GRO -Experiment Interface 


a. 

Height 

1830 kg 

b. 

Size 

2.25 m x 1.65 m diameter 

c. 

Power 

180 w including heaters 

d. 

Bit Rate 

6.8 kbps 


should be measurable on a scale of a few degrees. Features which sub- 
tend more than about 0.5° will be resolved above 1 GeV as extended 
sources. 

The diffuse radiation away from the galactic plane will be separable 
into galactic and extragalactic components on a scale of a few degrees. 
The extragalactic component will be studied for spatial variations in 
intensity and spectral distribution, providing important clues for 
determining the origin of this radiation. 

With respect to extragalactic discrete objects, the situation is not 
so easy to predict. The radio galaxy CEN A: two seyfert galaxies, NGC 
4151 and MCG8-11-11; and the quasar 3C273 have been seen in y-rays-. The 
high energy telescope on GRO will have about ten times lower threshold 
for the detection of extragalactic objects than COS-B which detected 
3C273. If the y-ray luminosity distribution is similar to that for 
either x-ray or visible photons, as many as 30 quasars might be expected 
to be seen. Even pessimistic estimates suggest that at least ten should 
be seen. Although even more difficult to predict, several normal 
galaxies should be detectable and for the nearest ones a study of their 
structure in y-rays will be possible. 

References 

Hughes, E.B. et al. (1980), IEEE Trans. Nuc. Sci. NS-27 , p. 364. 


spark chambers with Nal total 

energy counter 

12 to 30,000 MeV 

15X 30 MeV - 10,000 MeV 

6560 cm‘ 

2000 cnr above about 200 MeV 
5 x 10“®cm _2 s” 1 

1.6° at 100 MeV: 0.6* at 500 

MeV: 0.2° at 2 GeV (projected 

rms) 

5 arc min 

~. 5 sr: maximum opening angle 45° 
0.1 ms absolute 
2 years 
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A LARGE AREA EXPERIMENT TO 
DETERMINE COSMIC RAT ISOTOPIC ABUNDANCES 

B.G. Mauger* * V,K< Balasubrahmanyan, J.F. Ones, and R.F. StreitmatUi 
NASA/Goddard Space Flight Center, Greenbelt, MD 20771, U.S.A. 

W. Heinrich, Vi. Simon, H. 0. Tlttel 
Department of Physics, University of Siegen, 29 Slegen, West Germany 

1. ABSTRACT 

Observations of the cosmic ray isotopes of neon and magnesium have 
shown an overabundance of the neutron rich isotopes as compared with 
solar abundances* In order to determine the metalicity of the cosmic 
ray sources, the isotopic composition of these and other fluxes must be 
studied. A 1*2 m -sr instrument to study the Isotopic composition of 
the elements from oxygen through argon is being constructed as a joint 
project between the University of Siegen, West Germany and the Goddard 
Space Flight Center, U.S.A. This instrument uses two scintillators and 
two Cherenkov detectors (n ■ 1.33 and 1.4). Each of these radiators 
must have a uniform response to within 1 percent in order to achieve a 
mass resolution of .3 amu. To achieve this uniformity of response, each 
radiator will be frosted to eliminate specular reflections and will be 
used to determine the trajectory within ± 1 mm. A passive nitrocellu- 
lose range stack determines the range of the particles. These sheets of 
nitrocellulose are etched in an ammonia bath and the position of the 
tracks are read by computerized scanning. It is expected that this 
experiment will be able to determine the isotopic abundances in the 
energy range of 300 to 800 MeV/nucleon for elements from oxygen through 
argon. 

1 . INTRODUCTION 


Experiments to measure the isotopic composition of cosmic rays have 
shown an enrichment of neutron rich isotopes as compared with solar 
abundances. The cosmic ray isotope ratios, ZZ Ne/ Z0 Ne and ( Z5 Mg + 
z ®Mg)/ z Sig, exceed the solar abundance ratios by factors of 2.7 and 1.8, 
respectively (Fisher et al. , 1976; Greiner et al., 1979; Mewalt et al., 
1980). 

There are several processes that could be responsible for the , 
excess of these neutron rich isotopes. (For a discussion, see Woosley 
and Weaver, 1981.) The implications of these models also predict 
various amounts of neutron enrichment of other, less abundant, 
species. If measured, these Isotopic abundances could determine which 
processes dominate in cosmic ray sources. In order to measure the 
isotopic abundances of these species, a large area experiment with a 
mass resolution of better than .3 amu is necessary. 
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Such an experiment la under construction as a joint affort between 
the University of Slogan, West Germany and NASA's Goddard Space Flight 
Canter. This experiment, the Aluminum Isotopic Composition Experiment 
(ALICE), nses a Cherehkov-range technique to determine isotopic 
composition of elements from oxygen through argon. The experiment has a 
geometry factor of 1.2 m 2 st and uses two Cherenkov detectors to 
Increase its dynamic range. 4 

II. DETECTOR COMPLEMENT 


Figure 1 schematically shows the detector complement. Two 
scintillators, SI. and S2, form a coincidence telescope to define an 
event. These scintillators are 120 cm x 120. cm x 1 cm Polycast "PS— 10" 
scintillators. In order to provide a uniform response, the surfaces of 
these scintillators are frosted (to eliminate specular teflections) , 
housed in diffusive light collection boxes, and viewed by 16 12.5 cm 
photomultiplier tubes. The interior of the light collection boxes is 
painted with a barium sulfate paint with a .95 reflectivity. The light 
collection efficiency of these boxes is on the order of 50 percent. 

The trajectory of the cosmic ray is determined by the drift 
chambers, DCI and DCII. Prototype drift chambers have been tested in 
the heavy ion beam at Berkeley (Simon et al. , 1981). Ti.ese chambers 
should provide sufficient spatial resolution to trace the ccsiiic ray 
through the detector stack to i 1 mm. 

The two Chereknov detectors Cl and CII, are housed in light 
diffusion boxes identical to those containing the scintillators. 

Because Chereknov detectors are most sensitive near threshold, two 
Cherenkov detectors with different indices of refraction are used. Cl 
is a 2 cm thick piece of Pilot 425 with a refractive index of 1.33. CII 
is a silicon based rubber, GE RTV-615, with a refractive index of 1.4. 

Belov DCII is a passive range stack of cellulose nitrate, CN. The 
range stack has a depth of 15 g/cm 2 (11 cm) made of .03 g/cm 2 (250um) 
foils. After the flight, these foils will be etched and the resultant 
tracks transferred to a print by an ammonia scan method. These prints 
will then be digitized by a computerized scanning microscope. 

At the bottom of the range stack, there is an anticoincidence 
scintillator, ANTI, which is used to reject particles which are not 
stopped by the range stack. 

III. MASS RESOLUTION 


In order to achieve e mass resolution of better than .3 amu, it is 
necessary to determine the range and the energy per nucleon (from the 
Chereknov detectors) to within 1 percent. One source of range error is 
struggling, which cannot be reduced by experimental methods. Another 
source of range error occurs due to uncertainties in the amount of 
material between the velocity (Cherenkov) meesurement and the range 
stack. Special design consideration has been given to the construction 
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of tho light diffusion box** snd th« drift chambers to koop mss 
uncertain tlss below « fsr hundredths of s g/en 2 . 


Errors In tbs velocity measurement csn occur for a variety of 
reasons* To eliminate variations In the Chereknov response* the double- 
diffuse light gathering technique described shove Is used. 

i 

To remove any residual response nonunif oral ties, the detectors vlll 
be response-napped using relativistic partlclen before and during the 
flight. VI th a Mpped response profile* the unif oral ties of the detec- 
tors should be reduced to below 1 percent. Since the photo cathodes 
represent a saall fraction of the diffusion box's area, on the average* 
each ;pho ton suffers nany diffuse reflections. before being detected. 
Finally photon counting statistics Halt the accuracy of the Cherenkov 
detectors. Figure .2 shows the estimated aass resolution of the 
instrument. 
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Figure 1. The ALICE detector complement 
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A METHOD TO DETECT ULTRA HIGH ENERGY ELECTRONS USING 
EARTH MAGNETIC FIELD AS A RADIATOR 


S. A. Stephens* end V. K. Bala Subrahmanyan 
NASA/Godderd Spece Flight Center, Greenbelt, MD 20771 

1. INTRODUCTION 

Detection of high energy electrons above a few TeV, which lose 
energy rapidly through synchrotron and Inverse Compton processes, would 
give valuable Information on the distribution of sources and on the 
propagation of cosmic rays In the solar neighborhood* The influence of 
a strong nearby source would also be seen as a marked enhancement over 
the depleted electron spectrum. However, none of the existing experi- 
mental techniques could be used to measure the energy spectrum beyond a 
few TeV. Ve examine here the possibility of detecting electrons [1] 
above a few TeV by the photons emitted through synchrotron radiation In 
the Earth's magnetic field. As these photons are emitted in the 
direction of the electron within an extremely narrow cone, they would be 
Incident on a detector almost along a straight line. This pattern 
provides an unique signature for the identification of the electron. 

The mean energy of the emitted photons, being proportional to the square 
of the energy, provides reliability in the energy estimation even if a 
detector has low energy resolution. Thus, this technique would open up 
an Important window of energy for the study of primary cosmic ray 
electrons . 

2. SYNCHROTRON RADIATION OF ELECTRONS IN EARTH’S MAGNETIC FIELD 


The theory of synchrotron radiation is extensively discussed in the 
literature [2] . The spectral shape of photons emitted by an electron 
has a slope of -2/3 at low energies, 6lowly steepens at higher energies, 
and drops off exponetially beyond Eq ■ 6.56 x 10 2 B^E 2 MeV; B 
in Gauss and E Q in TeV. It can be shown that 70 percent of the photons 
emitted are spread over 2 decades in energy’. The mean energy of the 
photons is proportional to the square of the electron energy, and In 
order to cover m decades of electron energy, one needs an instrument 
capable of detecting photons over (m 2 + 2) decades in energy. 

Photons are emitted within a narrow cone in the direction of 
electrons . The radiation is polarized, and the angular distributions of 
photons along the two principal directions of polarization are 
different. We infer from the angular distribution of photons that 70 
percent of the photons are emitted within 0.65 x 10”' radian. This 
correspond to a spread of only 1 .3 ems on to a detector placed at a 
distance of 100 km, implying that the photons are distributed approxi- 
mately along a straight line. 

With a typical balloon altitude corresponding to a depth of- 
4 g cm” 2 of atmosphere over Palestine, Texas as the initial point on the 
trajectory of the electron, we have traced its trajectory in Earth's 
field using sixth degree simulation of the geomagnetic field [3,4]. 
Figure la shows an exaggerated view of the trajectory of an electron. 
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In this figure 0 Is ths arrival point , B's danota tha staps with stap 
langth As, and A's ara tha lntarcapts on tha dataetor plana by tanfants 
*1*1* *2 a 2**’* At P°* nt B's on tha trajectory. Tha vary llnas are 
indicative of tha photons emitted by tha electron. At aach stap, we 
hava ealculatad tha nav coordinates, B., tha radius of curvature, tha 
distance svapt by tha tangent At, and I ■ ^ At. Tha energy spectrum of 
photons > 20 kaV Is corrected for tha absorption of photons in the air 
between the detector and the point of emission* The total number of 
photons received over A l, the energy flov, the mean energy <E^>, and 
their angular distribution are determined. 

3. DETECTOR RESPONSE AND ENERGY ESTIMATION 


When an electron passes through the detector, the photons would be 
incident only on the convex side of the trajectory (Fig. la). This Is a 
special situation making the identification of the electron event very 
unique. Even if the electron does not intercept the detector, one can 
still identify the electron by the colinearity of the incident photons. 
The number of photons required to identify an event unambiguously 
is n >3. We have calculated the response area A. of a detector using 
this^criterion. We divide the detector into a number of strips (Fig. 
lb) of width Ab parallel to the arrival direction. If the electron 
arrives at 0, the detector strip would have 3 photons (billed circles) 
between 0 and D, and none between 0 and A. Therefore, the minimum 
length required by the detector is OD - i. As the electron moves 
towards A, the number of photons detected Increases, reaching a maximum 
at A, and then decreases to less than 3 photons beyond 0. Thus, the 
maximum distance up to which an electron can be detected is 0"D ■ l2. 

The response area of the strip is (i 9 - x \ b. 
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The calculated response of a circular detector la shown in Fit* 2 
for electrons arriving at a zenith angle of 45*. The detector la shown 
by thick circles. The thin curves trace the circumference of the 
response area for n„ ■ 5* The effect of increasing the value of ty to 
7 photons is shown by the dashed curve for North. The dramatic Increase 
of the response area with can be noticed from this figure at 20 TeV; 
this figure is scaled down by a factor of 2. The response area (a) 
varies with the arrival direction of electrons, (b) appears elongated 
along the diameter, indicating its dependence on the shape and size of 
the detector, (c) decreases as the value of n„ r is Increased, and (d) 
Increases with energy. 

Knowing B , the energy of an electron can be determined from the 
observed mean energy of photons. Since, the number of particles 
detected by an instrument in a direction (Z, V) is proportional to the 
collecting power CP(Z,$) in that direction, one can make use of the 
calculated CP to draw histograms of the estimated <E > . CP is the 
product of the projected A R perpendicular to the arrival direction, and 
the solid angle is that direction. Using this, we examine the viability 
of using an omnidirectional detector to study electrons. Fig. 3 shows 

histograms of <E^> using an omnidirectional 
detector of radius 150 ems for E 0 “ 10 TeV and 
n^ • 5 photons, with zenith angle cuts as 
Indicated. The peak in the distribution shifts 
to higher energies S6 one Includes larger zenith 
angles, because of the higher mean B . It can 
be shown that <E.^> distribution is Independent 
of n^ and the detector area. The distributions 
are asymmetric. The half widths at half 
maximum, in terms of the energy of the electron, 
are about 15 and 5 percent of the peak value 
towards lower and higher energies. For a 
typical experiment discussed here, we obtain a 
relation for the most probable energy as 
<E >m • 0.0225 E^’ 76 MeV. It may be noted that 
this relation is based on the assumption that 
the Instrument has 100 percent detection 
efficiency for photons > 20 keV. In actual 
practice this is not true, and one needs to 
Fig. 3 Histograms of Include both the detector efficiency and 

the energy response to obtain a similar relation. 

4. DETECTION OF ELECTRONS 



Ue have calculated the collecting power of an omnidirectional 
detector and examined its dependence on E Q , the area of the detector, 
and the selection criterion n y . It is found that as long as A is not 
very small and n y not very large, the collecting power of a detector is 
asymptotically proportional to Ae n^, where a - 1.43 and B - - 1.8. As 
the energy decreases, CP deviates ?r v Its linear dependence on E and 
decreases by an order of magnitude ^ a decrease of energy by a factor 
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of 2. This behavior provides the threshold detection for electrons, end 
thus, the integral flux can be determined accurately even by an 
instrument with limited energy resolution. 

For an ideal detector, capable of detecting photons > 20 kaV, we 
have calculated the expected rate of electron events using an electron 
spectrua of the type J(E ) * 500 electrons/ (a 2 sr s GeV). The 

event rate thus calculated Is shown in Fig. 4 as a function of the 
energy of electron for various values of n^ and A. ' Curves A, 8, and C 
In this figure correspond to circular detectors with radii 200 eat, 100 
cat, and 50 cat, respectively. The collection time In a typical balloon 
flight Is about 20 hours. Using a detector of area 10m , a total of 
about 13 electrons can be observed at energies above 5 TeV using n 
criterion of 4 photons; this number reduces to about 8.5 for n ■ 5 
photons. Fron a 14-day space shuttle flight ope can collect about 100 
events above 5 TeV and about 3 events above 10 1 * eV. As there Is no 
atmospheric absorption of X-rays, the detection threshold la decreased 
and the number of events would be larger. We have compared these event 
rates with that expected fron a proposed [5] shuttle experiment as shown 
by the curve TRIC; the expected number of events 1c only 0.5 electron 
above 5 TeV. 



Fig. 4. Expected event rate 
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